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Abstract: We analyze the maximum branching ratios for the Flavor Changing Neutral 
Current (FCNC) decays of the neutral Higgs bosons of the Minimal Supersymmetric Stan- 
dard Model (MSSM) into bottom quarks, h ^ bs {h = hP,H^,A^). We consistently 
correlate these decays with the radiative B- meson decays {h — > 57). A full-fledged com- 
bined numerical analysis is performed of these high-energy and low-energy FCNC decay 
modes in the MSSM parameter space. Our calculation shows that the available data on 
B{h 57) severely restricts the allowed values of B{h bs). While the latter could 
reach a few percent level in fine-tuned scenarios, the requirement of naturalness reduces 
these FCNC rates into the modest range B{h — > 5s) ~ 10~^ — 10"'^. We find that the 
bulk of the MSSM contribution to B(h bs) could originate from the strong supersym- 
metric sector. The maximum value of the FCNC rates obtained in this paper disagree 
significantly with recent (over-)estimates existing in the literature. Our results are still 
encouraging because they show that the FCNC modes h ^ bs can be competitive with 
other Higgs boson signatures and could play a helpful complementary role to identify the 
supersymmetric Higgs bosons, particularly the lightest CP-even state in the critical LHC 
mass region mf^o ~ 90 — 130 GeV. 
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1. Introduction 



Experimentally, processes involving Flavor Changing Neutral Current (FCNC) have been 
shown to have extremely low rates [1]. Theoretically, their rareness can be explained by 
the GIM mechanism [2], which is related to the unitarity of the mixing matrices between 
quarks. The Minimal Standard Model (SM) embeds naturally the GIM mechanism, due 
to the presence of only one Higgs doublet giving mass simultaneously to the down-type 
and the up-type quarks, and as a result no tree-level FCNCs interactions appear. FCNCs 
are radiatively induced, and are therefore automatically small. The addition of further 
Higgs doublets to the SM in the most general way introduces tree-level FCNC interactions, 
which would predict significant FCNC rates. However, by introducing an ad-hoc discrete 
symmetry these interactions are forbidden. This gives rise to two classes of Two-Higgs- 
Doublet Models (2HDM) which avoid FCNCs at the tree-level, known conventionally as 
type I and type II 2HDMs [3]. 

Supersymmetry (SUSY) [4-7], on the other hand, provides an appealing extension of 
the SM, which unifies the fermionic and bosonic degrees of freedom of the fundamental 
particles, and provides a natural solution to the hierarchy problem. The search for SUSY 
particles has been one of the main programs of the past experiments in high energy physics 
(LEP, SLD, Tevatron), and continues to play a central role in the present accelerator 
experiments (Tevatron II), and in the planning of future experimental facilities, like the 
LHC and the LC [8-10]. The Minimal Supersymmetric Standard Model (MSSM) is the 
simplest extension of the SM which includes SUSY, and for this reason its testing will be 
one of the most prominent aims of these powerful experiments. 

Complementary to direct searches, one can also look for effects of particles beyond 
the SM by studying their radiative effects. Much work along these lines has been made 
over the past two decades. FCNCs may play an important role here because they are 
essentially loop-induced. Hence SM and non-SM loops enter the FCNC observables at the 
same order of perturbation theory, and new physics competes on the same footing with 
SM physics to generate a non-vanishing value for these rare processes. It may well be that 
the non-SM effects are dominant and become manifest. Conversely, it may happen that 
they become highly constrained. The power of FCNC observables can be gauged e.g. by 
the implications of the bottom-quark rare decay b sj: the experimentally measured 
allowed range B{b — 57) = (3.3 it 0.4) x 10~^ [1, 11-16] may impose tight constraints on 
extensions of the SM. For example, it implies a lower bound on the charged Higgs boson 
mass mH± > 350 GeV in general type II 2HDMs [17-20]. 

The most general MSSM includes tree-level FCNCs among the extra predicted parti- 
cles, which induce one-loop FCNC interactions among the SM particles. Given the observed 
smallness of these interactions, tree-level SUSY FCNCs are usually avoided by including 
one of the two following assumptions: either the SUSY particle masses are very large, and 
their radiative effects are suppressed by the large SUSY mass scale; or the soft SUSY- 
breaking squark mass matrices are aligned with the SM quark mass matrix, so that both 
mass matrices are simultaneously diagonal. However, if one looks closely, one soon realizes 
that the MSSM does not only include the possibility of tree-level FCNCs, but it actually re- 
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quires their existence [21]. Indeed, the requirement of SU{2)l gauge invariance means that 
the up-left-squark mass matrix can not be simultaneously diagonal to the down-left-squark 
mass matrix, and therefore these two matrices can not be simultaneously diagonal with the 
up-quark and the down-quark mass matrices, that is, unless both of them are proportional 
to the identity matrix. But even then we could not take such possibility too seriously, for 
the radiative corrections would produce non-zero elements in the non-diagonal part of the 
mass matrix. All in all, we naturally expect tree-level FCNC interactions mediated by the 
SUSY partners of the SM particles. The potentially largest FCNC interactions are those 
originating from the strong super symmetric (SUSY-QCD) sector of the model (viz. those 
interactions involving the squark-quark-gluino couplings), and in this paper we mainly 
concentrate on them. These couplings induce FCNC loop effects on more conventional 
fermion-fermion interactions, like e.g. the gauge boson-quark vertices Vqq' . 

Of course, low energy meson physics puts tight constraints on the possible value of 
the FCNC couplings, especially for the first and second generation squarks which are 
sensitive to the data on — and — [22,23]. The third generation system is, 
in principle, much loosely constrained, since present data on — mixing still leaves 
a large room for FCNCs, and the most stringent constraints are given by the B{h — > 57) 
measurement [11-16]. Therefore the relevant FCNC gluino coupling 822, [22,23] (see Section 
HI) is not severely bound at present. The lack of tight FCNC constraints in the top-bottom 
quark doublet enables the aforementioned lower bound on the charged Higgs boson mass in 
the MSSM to be easily avoided, to wit: by arranging that the SUSY-electroweak (SUSY- 
EW) contribution to B{h ^ sj) from the top-squark/chargino loops screens partially the 
charged Higgs boson contribution. This situation can be naturally fulfilled if the higgsino 
mass parameter (/x) and the soft SUSY-breaking top-squark trilinear coupling (At) satisfy 
the relation fiAt <0 [24-28]. 

The FCNC gluino interactions also induce large contributions to i?(6 — > 57). It should 
however be noted that the leading contributions to the Vqq' FCNC interactions from the 
third quark generation correspond to a double insertion term, in which the squarks propa- 
gating in the loop suffer a double mutation: a flavor conversion and a chirality transition. 
This fact has been demonstrated in the S(6 — > 57) observable itself [29], as well as in 
the FCNC rare decay width T{t eg) [30,31]. As a consequence, the loose limits on the 
third generation FCNC interactions derived under the assumption that the leading terms 
contributing to 6 ^ 57 correspond to the single particle insertion approximation [22, 23] 
are not valid, and more complex expressions must be taken into account [32]. 

Concerning the FCNC interactions of Higgs bosons with third generation quarks, it 
was demonstrated long ago [30, 31] that the leading term corresponds to a single parti- 
cle insertion approximation, which produces a flavor change in the internal squark loop 
propagator, since in this case the chirality change can already take place at the squark- 
squark-Higgs boson interaction vertex. Adding this to the fact that the Higgs bosons (in 
contrast to gauge bosons) have a privileged coupling to third generation quarks, one might 
expect that the FCNC interactions of the type quark-quark-Higgs bosons in the MSSM 
become highly strengthened with respect to the SM prediction. This was already proven 
in the rare decay channels T(t — > ch) [30,31] {h being any of the neutral Higgs bosons of 
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the MSSM h = h^,H°,A°), where the maximum rate of the SUSY-QCD induced bran- 
ching ratio was found to be BR{t — > ch) ~ 10~^, eight orders of magnitude above the SM 
expectations BR{t — > cH^^^) ~ 10^^^. Similar enhancement factors have been found in 
the top-quark-Higgs boson interactions in other extensions of the SM [31,33,34]. 

From the experience of the previous calculations with the top quark, we expect similar 
enhancements in the FCNC interactions of the MSSM Higgs bosons with the bottom quark. 
Indeed, the purpose of this paper is to quantify, in a rehable way, the MSSM expectations 
on the FCNC Higgs boson decay modes 

h^hs, h^bs {h = h^,H^,A^). (1.1) 

There are other FCNC decay modes involving light quarks. However, only these bottom 
quark channels are relevant, as the remaining FCNC decays into light quarks are negligible 
in the MSSM. Moreover, the FCNC decays of Higgs bosons into bottom quarks are specially 
interesting as they can provide an invaluable tool to discriminate among different extended 
Higgs boson scenarios in the difficult LHC range 90 < mh < 130 GeV [8,9]. 

In this paper we present what we believe is the first realistic estimate of the SUSY- 
QCD contributions to the FCNC branching ratios of the MSSM Higgs bosons into bottom 
quark. Specifically, we compute 

for the three Higgs bosons of the MSSM, h = h^,H°,A°, where T{h X) is the - 
consistently computed - total width in each case. The maximization process of the above 
branching ratios in the MSSM parameter space is performed on the basis of a simultaneous 
analysis of the relevant partial decay widths and of the branching ratio of the low-energy 
FCNC process 6 — > 57, whose value is severely restricted by experiment [11-16]. It turns 
out that the maximum FCNC rates that we find disagree quite significantly with some 
simplified estimates that have recently appeared in the literature [35]. According to these 
authors the FCNC decay rate of some MSSM Higgs bosons into bottom quarks can reach 
the level B(h qq') ^ 25%. We find this value untenable, even more given the fact that 
in Ref. [35] no attempt is made to verify the restrictions of the parameter space imposed 
by the low energy data on i?(6 — > 37).^ 

The structure of the paper is as follows. In Section ^ we estimate the expected bran- 
ching ratios and describe the structure of Eq. ( |1.2| ) in the MSSM in more detail; in Section 
m we present the numerical analysis, and in Section |^ we deliver our conclusions. 

2. Partial widths and branching ratios 



First of all let us estimate the branching ratio ( [1.2| ) in the SM. It is not necessary to perform 
a detailed calculation to suspect that it is rather small. Using dimensional analysis, power 



^See also Ref. [36] for a combined analysis of flavor-violating and CP-violating MSSM couplings. 
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counting, CKM matrix elements and dynamical features we expect that the maximum 
branching ratio of the SM Higgs boson into bottom quark is at most of order ^ 

BR{H^^' ^h-s)^(^^^\wGF(^^<lQ-' {iimH<2Mw). (2.1) 

Here Gp is Fermi's constant and aw = (7^/47r, g being the SU{2)l weak gauge coupling. 
This result should hold for a SM Higgs boson mass ran < 2MW, and in particular in the 
critical LHC region mn — 90 — 130 GeV. We have approximated the loop form factor by 
just a constant prefactor, and the numerical value is an approximate upper bound assuming 
that we approach mn = 2M]y from below. We have taken Vts = 0.04 and nib = 5 GeV for 
this estimate. In spite of the crudeness of the estimate, direct evaluation with programs 
FeynArts, FormCalc and LoopTools [37-39] confirms the order of magnitude ( |2.lD ^. On 
the other hand, if uih > 2M\y, more specifically if uih > rrit, it is easy to see that ( p.l| ) 
will be suppressed by an additional factor of ml/mjj, because the vector boson Higgs decay 
modes — > W~{Z Z) will be kinematically available and become dominant. But 
at the same time the ratio m^/m^ is replaced by m^/m^. Hence, 

BR{H^^^ ^ 6 S-) ~ (^11^) ' aw Gf < IQ-^o ( for mn > nit) . (2.2) 

In the numerical evaluation we assumed a mass range where the ratio mn/'mt > 1 is of 
order one as this provides an upper bound. In both cases (2.1) and ( |2.2D the branching 
ratios into bottom quark are much larger than the Higgs boson FCNC branching ratio into 
top quark in the SM [34]. However, even in the case ( p.lD it is still too small to have 
a chance for detection in the LHC. It is clear that unless new physics comes to play the 
process course H^^^ — > 6s) will remain virtually invisible. Nonetheless 

the result ( |2.1| ) is not too far from being potentially detectable, and one might hope that 
it should not be too difficult for the new physics to boost it up to the observable level. 

Consider how to estimate the potentially augmented rates for the MSSM processes 
( |1.1| ), if only within a similarly crude approximation as above. Because of the strong 
FCNC gluino couplings mentioned in Section |l] and the tan /3-enhancement inherent to 
the MSSM Yukawa couplings (see Ref. [30,31] for details), we may expect several orders 
of magnitude increase of the branching ratios ( |1.2D as compared to the previous SM re- 
sult. A naive approach might however go too far. For instance, one could look at the 
general structure of the couplings and venture an enhancement factor typically of order 
(as/aw)'^ tan^ /? |523/Vts|^, which for 623 < 1 and tan/3 > 30 could easily rocket the SM 
result some 5 — 6 orders of magnitude higher, bringing perhaps one of the MSSM rates 
( |1.2| ) to the "scandalous" level of 10% or more. But of course only a more elaborated 
calculation, assisted by a judicious consideration of the various experimental restrictions. 



^See Ref. [34] for details on similar estimates, like that of BR{H^'^' ~* tc). 

^To our knowledge, the first evaluation of the SM branching ratio was performed in [40], however no 
detailed analysis of it exists in the literature. It is natural to clarify this value before jumping to evaluate 
the possible non-SM contributions. For the details, see Ref. [41]. 
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can provide a reliable result. As we shall see, a thorough analysis generally disproves the 
latter overestimate. 

The detailed computation of the SUSY-QCD one-loop partial decay widths T{h — > q q') 
in (|l3) within the MSSM follows closely that of r(t ch) (see Ref. [30,31]). The rather 
cumbersome analytical expressions will not be listed here as they are an straightforward 
adaptation of those presented in the aforementioned references. However, there are a few 
subtleties that need to be pointed out. One of them is related to the calculation of the 
total widths T{h X) for the three Higgs bosons h = h?, , in the MSSM. As long as 
r(/i — > g q') in the numerator of Eq. ( |l.2D is computed at leading order, the denominator has 
to be computed also at leading order, otherwise an artificial enhancement of B{h q q') can 
be generated. For example, including the next-to-leading (NLO) order QCD corrections to 
T{h bb) reduces the decay width by a significant amount [42-46]. Then, to be consistent, 
the NLO (two-loop) contributions to T{h qq') should also be included. Similarly, 
the one-loop SUSY-QCD corrections to T{h — > bb) can be very large and negative [47], 
which would enhance B{h — > qq'). At the same time these corrections also contribute 
to r(/i qq'), such that contributions to the numerator and denominator of Eq. ( |1.2| ) 
compensate (at least partially) each other. Therefore the same order of perturbation theory 
must be used in both partial decay widths entering the observable B{h — > qq') to obtain 
a consistent result. By the same token, using running masses in the numerator of ( |l.2| ) 
is mandatory, if they are used in the denominator. Last, but not least, consistency with 
the experimental bounds on related observables should also be taken into account. In 
this respect an essential role is played by the constraints on the FCNC couplings from 
the measured value of B{b 57). They must be included in this kind of analysis, if we 
aim at a realistic estimate of the maximal rates expected for the FCNC processes ( |l.]j ) 
in the MSSM. In our calculation we have used the full one-loop MSSM contributions to 
i?(6 — > 57) as given in [48] ^. 

Let us now summarize the conditions under which we have performed the computation 
and the approximations and assumptions made in the present analysis: 

• We include the full one-loop SUSY-QCD contributions to the partial decay widths 
T{h^qq') in(|L2D. 

• We assume that FCNC mixing terms appear only in the LH-chiral sector of the squark 
mixing matrix. This is the most natural assumption, and, moreover, it was proven in 
Ref. [30, 31] that the presence of FCNC terms in the RH-chiral sector enhances the 
partial widths by a factor two at most - not an order of magnitude. 

• The Higgs bosons total decay widths r(/i — > X) are computed at leading order, 
including all the relevant channels: T(h — > //, ZZ, W~^W~ , gg). The off-shell decays 
T{h Z Z* ,W^W^*) have also been included. The one-loop decay rate T{h gg) 
has been taken from [49] and the off-shell decay partial widths have been computed 
explicitly and found perfect agreement with the old literature on the subject [50]. 

*Ref. [48] contains a partial two-loop computation of B(b ~* 57), but only the one-loop contributions 
have been used for the present work. 
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We have verified that some of the aforementioned higher order decays are essential 
to consistently compute the total decay width ofT{h^ — > X) in certain regions of the 
parameter space where the maximization procedure probes domains in which some 
(usually leading) two-body processes become greatly diminished. We have checked 
that our implementation of the various Higgs boson decay rates is consistent with the 
results of HDECAY [51]. However, care must be exercised if using the full-fledged result 
from HDECAY. For example, it would be inconsistent, and numerically significant, to 
compute the total widths r(/i — > X) with this program and at the same time to 
compute the SUSY-QCD one-loop partial widths T{h — > qq') without including the 
leading conventional QCD effects through e.g. the running quark masses. 

• The Higgs sector parameters (masses and CP-even mixing angle q) have been treated 
using the leading rrit and tan f3 approximation to the one-loop result [52-55] . For 
comparison, we also perform the analysis using the tree-level approximation. 

• We include the constraints on the MSSM parameter space from B{b ^ S7). We adopt 
B{b 57) = (2.1 — 4.5) X 10~^ as the experimentally allowed range within three 
standard deviations [1]. Only the SUSY-QCD contributions induced from tree-level 
FCNCs are considered in the present work. 

Running quark masses {mq{Q)) and strong coupling constants {as{Q)) are used throughout. 
More details are given below, as necessary. 



3. Full one-loop SUSY-QCD calculation: Numerical analysis 

Given the setup described in Section ^, we have performed a systematic scan of the MSSM 
parameter space with the following restrictions: 

523 < 10-0 09 ~ 0.81 

= -1500- •• 1500 GeV 
12 = -1000- •• 1000 GeV ^ ' -* 

= 150 • • • 1000 GeV 

and the following fixed parameters: 

tan/3 = 50 

At = -300 GeV . 

Here m^^ ^ are the left-chiral and right-chiral bottom-squark soft-SUSY-breaking mass 
parameters, and rriq is a common mass for the strange- and down-squark left- and right- 
chiral soft SUSY-breaking mass parameters. Following the same notation as in [30], the 
parameter 623 represents the mixing between the second and third generation squarks. Let 
us recall its definition: 

523^^^, (3.3) 
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Figure 1: Maximum SUSY-QCD contributions to B{h — > qq'), Eq. ( |l.2| ), as a function of a) m^o 
and b) ,523 for m^o = 200 GeV. 
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Table 1: Maximum values of B{h — > qq') and corresponding SUSY parameters for m^o = 200 GeV. 



m? _ being the non-diagonal term in the squark mass matrix squared mixing the second 
and third generation left-chiral squarks. The parameter is a fundamental parameter in 
our analysis as it determines the strength of the tree-level FCNC interactions induced by 
the supersymmetric strong interactions, which are then transferred to the loop diagrams 
of the Higgs boson FCNC decays 

The result of the scan is depicted in Fig. ||. To be specific: Fig. ||a shows the maximum 
value B™^^{h qq') of the FCNC decay rate ( ^^ ) under study as a function of m^o; 
Fig. displays B'^^^{h — > qq') as a function of the mixing parameter for rrij^o = 
200 GeV. Looking at Fig. |^ three facts strike the eye immediately : i) the maximum is 
huge (13%!) for a FCNC rate, actually it is as big as initially guessed from the rough 
estimates made in Section |2|; ii) very large values of 623 are allowed; iii) the maximum rate 
is independent of the pseudo-scalar Higgs boson mass m^o. We will now analyze facts 
ii) and iii) in turn, and will establish their incidence on fact i). For further reference, in 
Table [l| we show the numerical values of B'^''^^{h qq') together with the parameters 
which maximize the rates for m^o = 200 GeV. 
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Figure 2: B{h qq') and B{b sj) as a function of 623 for tlie parameters that maximize 
B{h^ ^ qq') in Table |l|. The shaded region is excluded experimentally. 
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Figure 3: Maximum value of the SUSY-QCD contributions to B{h qq') as a function of a) 
rriAO and b) ^23 for 171^0 = 200 GeV, for the scenario excluding the window regions. 



One would expect that a large value of 623 should induce a large gluino contribution 
to B(b — > 57). In fact it does! However our automatic scanning process picks up the 
corners of parameter space where the gluino contribution alone is much larger than the 
SM contribution, but opposite in sign, such that both contributions destroy themselves 
partially leaving a result in accordance with the experimental constraints. We examine 
this behaviour in Fig. |2|, where we show the values of B(h ^ q q') together with B{h ^ 57) 
as a function of ^23 for the parameters which maximize the FCNC rate of the lightest CP- 
even state /i° in Table |. We see that, for small values of 523 5 the gluino contribution to 

— > 57) is small, and the total B{b ^ 57) prediction is close to the SM expectation. In 
contrast, as 823 steadily grows, B{h — > S7) decreases fast (meaning a dramatic cancellation 
between the two contributions) until reaching a point where B[h sj) = 0. From there on 
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Table 2: Maximum values of B{h q q') and corresponding SUSY parameters for m^o — 200 GeV 
excluding the window region. 

it starts to grow with a large slope, and in its race eventually crosses the allowed — > 57) 
region. The crossing is very fast, and so rather ephemeral in the 623 variable, and it leads 
to the appearance of a narrow allowed window at large ^23 values, see Fig. |2|a. We would 
regard the choice of this window as a fine-tuning of parameters, hence unnatural. For 
this reason we reexamine the B{h — > qq') ratio by performing a new scan of the MSSM 
parameter space in which we exclude the fine-tuned (or window) region. The result for 
m^o = 200 GeV can be seen in Table |2| and Fig. ^. This time we see that the maximum 
values of B{h qq') are obtained for much lower values of 623, and the maximum rates 
have decreased more than one order of magnitude with respect to Table |l|, reaching the 
level of few per mil. These FCNC rates can still be regarded as fantastically large. Had 
we included the SUSY-EW contributions to B{b ^ 57), further cancellations might have 
occurred between the SUSY-EW and the SUSY-QCD amplitudes. Even more: since each 
contribution depends on a separate set of parameters, one would be able to find a set 
of parameters in the SUSY-EW sector which creates an amplitude that compensates the 
SUSY-QCD contributions for almost any point of the SUSY-QCD parameter space [41]. 
But of course this would be only at the price of performing some fine tuning, which is not 
the approach we want to follow here. 

On the other hand further contributions to B{b — > 57) might exist. In the most general 
MSSM, flavor-changing interactions for the right-chiral squarks {523rr), and mixing left- 
and right-chiral squarks (523L_r) can be introduced. The latter can produce significant con- 
tributions to B{b — > 57), changing the allowed parameter space. The introduction of 523lr 
can produce two possible outcomes: First, in certain regions of the parameter space, the 
contributions of 523LR and ^23 are of the same sign, enhancing each other. In this situation, 
the maximum allowed value of ^23 is obtained for 623LR = 0. Second, in other regions of the 
parameter space the two contributions would compensate each other, producing an overall 
value of B{b 57) in accordance with experimental constraints, even though each contri- 
bution would be much larger. Again, we would regard these compensations as unnatural, 
and would discard that region of the parameter space. In the following we will require that 
the SUSY-QCD contributions induced by ^23 do not compensate the SM ones to give an 
acceptable value of B{b — > 57); this is equivalent to the condition that the SUSY-QCD 
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Figure 4: B(h qq') and T{h^ X) (in GeV) as a function of /i for a) one-loop a angle; b) 
tree-level a angle, and for the parameters that maximize B{h" qq') in Table 0. The H° and A° 
curves coincide. The B{b ^ 57) constraint is not shown. 



amplitude represents a small contribution to the total B{h ^ 57) value, and is therefore 
independent of the inclusion of the other contributions (SUSY-EW, 52i,lr)-^ 

We turn now our view to the second fact, namely the independence of the maximum 
rates with respect to m^o . We will show that it also plays a central role as to the enhance- 
ment of B[h — > qq'). Actually, a good hint is given by the small values of the lightest 
Higgs boson decay width in Tables || and |2|, r(/i° ^ X) ~ 2 x 10^^ GeV. The maximiza- 
tion process of B{hP qq') does not only find the parameters for which T{hP qq') is 
maximum, but also the parameters for which T{h^ X) is minimum. Specifically, since 
T{h? bb) is the dominant decay decay channel of for large tan/3, the maximum of 
B{h^ qq') is produced in the parameter range of the so-called small a^jf scenario [56], 
that is, a parameter range where the radiative corrections make the CP-even Higgs boson 
mixing angle a vanish (or very small), such that the leading partial decay width T{h^ — > bb) 
is strongly suppressed. The consequences of this scenario have been extensively studied 
in Ref. [57]. As advertised in Section ||, the possibility that the maximization process 
explores these regions of the parameter space is the reason why the leading higher order 
decay channels, and also the leading three-body decay modes have to be taken into account 
in the computation of the total width. 

In Fig. ^ we plot the value of the various branching ratios B{h ^ q q') and of the total 
width of the lightest CP-even Higgs boson, T{h? X), as a function of the higgsino mass 
parameter ^, the rest of the parameters being those of the third column of Table |^, i.e. 
the ones that maximize the branching ratio B{hP qq')- Fig. ^a shows that T{h^ X) 
has a deep minimum in the range of corresponding to the maximum of B{h^ — > qq'), 
which reaches the level of a few percent. If, instead of using the radiatively corrected a 

^The analysis of Ref. [36] follows the opposite approach, that is: to find the fine-tuning conditions 
imposed by low energy data that allow for the largest possible value of the FCNC parameters. 
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4.50 X 10"^ 
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Table 3: Maximum values oi B{h ^ qq') and corresponding SUSY parameters for rn^io = 200 GeV, 
using the tree-level expressions for the Higgs sector, and excluding the window region. 

value we use the tree-level expression, we obtain the result shown in Fig. ^3. Here the 
total decay width of the Higgs boson is independent of /x, and B{h qq') does not show 
any peak. Actually in this case the branching ratio for /i" becomes smaller than that of 

and for all /j,. The maximization procedure in Fig. |^ selects for each value of m^o 
the MSSM parameters corresponding to the small Oe// scenario for that specific value of 
m^o. Of course, this discussion regarding the channels for large values of m^o has a 
correspondence with the channel for low values^ of my^o. 

As indicated in Section ^, we have used a one-loop approximation for the Higgs sec- 
tor [52-55], instead of the more sophisticated complete two- loop result present in the lite- 
rature [58,59]. However, we should stress that the exact MSSM parameters at which the 
small cxeff scenario is realized are not important for the sake of the present analysis. All 
that matters is that some portion of the parameter space exists, for which T{h^ bb) is 
strongly suppressed, but r{h^ — > qq') is not. 

To compare the maximum value of B(h^ qq') obtained with and without the small 
Oieff scenario, we have performed the maximization procedure using the tree- level expres- 
sions for the Higgs sector parameters. The result is shown in Table |3| and Fig. ^. In this 
case B{h^ — > qq') is reduced by a sizeable factor of > 20 with respect to Table |2|, whereby 
the rate descends about an order of magnitude below that of the /A^ channels which 
remain basically unchanged. Notice also that T{h^ X) is larger than in previous tables. 
In spite of the reduction, achieving a FCNC ratio B(h^ ^ (7 g') ~ 1.3 x lO""* is a remarkable 
result, three orders of magnitude larger than the maximum SM rate ( |2.1| ), and only one 
order of magnitude below the rare decay B{h^ — > 77) ~ 10~^. Also worth noticing in 
Fig. ^3 (and Fig. ^3) is the fact that B^^^{h — > qq') is essentially flat in 523 in the upper 
range down to 623 ~ 10"°-^ ~ 0.16. The reason lies in the correlation between B{h ^ qq') 
and B{b — > 57). In order to comply with the (non-fine-tuned) value of B{b 57) for large 
823, the absolute value of the fi parameter must be small. When ^23 decreases, can grow 
to larger values, leaving the overall maximum rates B'^^^{h qq') effectively unchanged 
(see Eq. (|3.5|) below). 

^Large or low values here means m^o > n^JJlf or m^o < m'l^o^, i.e. above or below the maximum 
possible value for the mass of the lightest Higgs boson /i", respectively. 
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Figure 5: Maximum value of the SUSY-QCD contributions to B{h — > qq') as a function of a) 
m^o and b) 523, for m^o = 200 GeV and for the scenario excluding the window region and using 
the tree-level expressions for the Higgs sector parameters. 



The maximization process selects a squark mass scale in the vicinity of the maximum 
values used in the scanning procedure. We should point out, however, that the same order 
of magnitude for B{h — > qq') could be obtained with a much lower squark mass scale. In 
this case the lighter squark masses induce a much larger B(6 — > 57) value, and S23 is much 
more constrained. For example, if we perform a scan in the parameter space but 
fixing the squark mass scale to be niq < 500 GeV, we obtain the following values for the 
maximal branching ratios for niy^o = 200 GeV: 

^max^^o ^qq') = 1.4 X 10-^ , S"^^^(FV^° -^qq')= 9.2 x 10"^ , (3.4) 

with 823 ~ lO"'''^, /i ~ —110 GeV, and we have limited ourselves to the scenario avoiding 
the window regions and using the tree-level expression for the Higgs sector parameters. 
These numbers have to be compared with Table ^. 

The reason behind this scale independence admits an explanation in terms of an effec- 
tive Lagrangian approach [41], in which one can estimate the leading effective coupling to 
behave approximately as: 



y/2Mw cos (3 37r M|^c;y 



sin(/?-a) {H^) 

cos(/?-a) (/jO) . (3.5) 
sin(2/3) (^0) 



Aside from ensuring (at least) a partial SUSY scale independence of the leading terms, 
this expression also shows that B[h — > qq') has a weak dependence on the soft-SUSY- 
breaking trilinear coupling Ai,. The observed situation is similar to the flavor-conserving 
hhh interactions, where the cancellation of the A^, terms at leading order has been recently 
proven [60]. It also shows that the leading non-decoupling SUSY contributions to T{h^ — > 
q q') eventually fade out as the decoupling limit of the Higgs sector is approached: cos(/? — 
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Figure 6: B{h qq') and B{b — > 57) as a function of rrig for the parameters that maximize 
B{h'^ — > bs) exchiding the window region (see third column of Table The shaded region is 
excluded experimentally. 
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Figure 7: As in Fig. ^ but including the window region. The remaining parameters are fixed as 
in the third column of Table |l]. 



a) 0. We have found (using the tree-level expression for a) that the non-leading (SUSY- 
decoupling) contributions to T{h^ qq') dominate for m^o > 450 GeV, inducing a value 
pmax^^jO ^ gg/) ^ 12 X 10-5, with 623 ~ W'^ , fi ~ 1000 GeV. Full details on the effective 
Lagrangian approach, and its application to further refine these calculations, will be given 
in a forthcoming publication [41] . 

We further investigate the role of the scale of SUSY masses, and the fine-tuning be- 
haviour in Figs. ^ and 0. In these figures we give up the equality rrig = ruq ( |3.2[) , the squark 
masses are fixed at the values stated in Tables |2| and || respectively. Fig. ^ shows the values 
of B{h — > qq') for the three Higgs decays and of B{h — > 57) as a function of the gluino 
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Table 4: Maximum values of B{h qq') and corresponding r{h X) for the different scenarios 
studied in this work. 



mass for the parameters that maximize B{h^ q q') when the window regions are excluded 
(third column of Table Q). Here we see that, while the gluino contribution to — > 57) 
decouples fast as a function of rUg, its contribution to B{h qq') is fairly sustained. 
Indeed, between rUg = 1 TeV and rrig = 5 TeV B{h^ — > qq') decreases only by a factor 
~ 1/4, while the gluino contribution to B{b — > 57) becomes negligible at nig = 5 TeV and 
we recover the SM prediction. As a consequence, the maximum rates B{h — > qq') that we 
have found are robust, in the sense that further theoretical refinements and experimental 
results that change the allowed range of B{b — > sj) can easily be compensated for by a 
slight increase of the gluino mass (mg), which would leave the prediction for B{h qq') 
essentially unchanged. We note in Fig. ^ the corresponding behaviour of B{h^ qq') and 
B{h ^ 57) in the presence of fine-tuning, i.e. as in Table |I[ In contrast to the previous 
case, here we observe the presence of two tiny windows in the regions rrig = 25 — 75 GeV 
and rrig = 950 — 1125 GeV. In the middle region rrig = 75 — 950 GeV, B{b ^ 57) is one 
order of magnitude larger than the allowed experimental range, and in the region above 
rrig = 1125 GeV it only enters the allowed region for nig > 4500 GeV. In this region 
B{h^ qq') is still large, but at the price of having a gluino five times heavier than the 
rest of the SUSY spectrum. This is another manifestation of the large fine-tuning that 
governs this region of the parameter space. 

Up to this point we have used the high tan/3 value quoted in Eq. ( |3.2| ). But we have 
also looked at the impact of varying tan/3 on B"^^^{h qq')- Since the latest LEP data 
restricts tan/? > 2.5, we have used a moderate value of tan/3 = 5. Note that, at low 
tan/3, the small Oe// scenario does not arise. As a consequence similar results are obtained 
using either the tree-level or one-loop expressions for the Higgs sector parameters. We 
find that the three branching ratios B"^^^{h — > qq') at tan/3 = 5 stay in the same order 
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of magnitude as in the scenarios with tan /3 = 50 (default case) with the tree-level Higgs 
sector and no-window (Cf. Table P). 

4. Remarks and conclusions 

The main numbers of our analysis are put in a nutshell in Table ^, where we show the 
results presented previously, together with some other scenarios and the low tan (3 case. 
The computed maximum values of B{h — > qq') must not be taken as exact numbers in 
practice, but order of magnitude results. The implications that can be derived from Table ^ 
can be synthesized as follows: 

1. The SUSY-QCD contributions can enhance the maximum expectation for the FCNC 
decay rates B{h qq') enormously. This is seen by comparing the results of Table § 
with the maximum value of B{H^^^ bs) considered in Eq. (|2.1| ). The optimized 
MSSM branching ratios are at the very least 3 orders of magnitude bigger than the 
SM result. 

2. If no special circumstances apply, that is, if no fine-tuning occurs between the parame- 
ters contributing to B{b — > 57) in the MSSM, and if r{h^ — > bb) is not suppressed, the 
maximum rates are qq') 1.3 x 10"^, ^ gg') ~ 9 x 10"^. 
This corresponds to the tree-Higgs / no-window scenario in Table |. 

3. If, however, T(h^ bb) is suppressed by the radiative corrections to the CP-even 
mixing angle a, then B{h^ qq') can be an order of magnitude larger: B™^^{h? — > 
qq') ~ 3 X 10^'^. This corresponds to the small Oe// scenario, and is indicated by 
small-ag/// Tio-window in Table ^. The FCNC branching ratio that we find for h? 
in this case should be considered as the largest possible one within the conditions of 
naturalness (no fine-tuning). 

4. On the other hand, if fine-tuning between the gluino and the SM contributions to 
B{b sj) is allowed, but the small-aejj scenario is not realized, then B™^^{hP — > 
qq') grows one order of magnitude up to B™^^{h^ qq') ~ 4 x 10~^, whereas 
B^^^{H^ /A^ qq') ~ 3 x 10~^. This corresponds to the case labelled tree- 
Higgs /window in Table ^. 

5. When both special conditions take place simultaneously, viz. fine-tuning in B{b 
s'j) (triggered by a very special choice of the ^23 parameter in a narrow window 
range) and small Ueff scenario (independent of assumptions on 623), we reach an 
over-optimistic situation where B™^^{hP — > qq') could reach the ~ 10% level. This 
is the case referred to as small-Og/ / / window in Table ^. 

6. If tan/3 is low/moderate, then B™^^{h — > qq') lie in the lower range ~ 10~^, which 
can grow an order or mag nitude for 5™''''(iJ°/y4° qq') in fine-tuned scenarios (last 
three rows in Table ^). 
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Although the large FCNC rates mentioned in points ^ and |5| above seem to offer 
a rather tempting perspective, we will not elaborate on them any further since in our 
opinion the fine-tuning requirement inherent in them is too contrived. On the other hand, 
points ^ and ^ offer a moderate, but certainly much more realistic scenario, which in no 
way frustrates our hopes to potentially detect the FCNC Higgs boson decays (1.1). Indeed, 



in the case described in point ^, B{h —>■ qq') can be at most of order 10~^. But this is still 
a fairly respectable FCNC branching ratio (comparable to that of 6 — > 57) and it may lead 
to a large number of events at a high luminosity collider [41]^. Moreover, if T{h X) 
becomes suppressed (e.g. by realizing the small Oeff scenario, point ^) then B{h^ —>■ qq') 
can be enhanced by an additional order of magnitude. 

Our analysis correlates the values of B(h qq') with that of B(b — > 57), taking into 
account only the SUSY-QCD contributions due to flavour mixing parameters among the 
left-chiral squarks. The presence of several other competing contributions to B(b — > 57) 
alters the borders of the allowed parameter space: 

• For the fine-tuned scenarios, the presence and position of the allowed window regions 
in the parameter space depends significantly on all the contributions, and therefore 
also does the maximum value of B{h — > qq'). Outside the window regions, the 
computed value B{b 57) can only be made consistent with the experimental 
range, by means of a large splitting between the squark and gluino masses. 

• For the non- fine-tuned scenarios, the inclusion of further contributions to B{b ^ 57) 
also alters the allowed parameter space, but the condition of non-fine-tuning ensures 
precisely that the change in the allowed range of 823 is smooth, and the corresponding 
change in B^^^{h — > qq') is not dramatic. 

Of course, the question immediately arises on what will happen if the data from present 
-B-meson factories further constrains the 823 parameter. In that case, we should take into 
account the (charged-current induced) SUSY-EW contributions to B(h qq'), which will 
be presented in Ref. [41] (see also [61]). However, we can advance that the SUSY-EW 
effects on B^^^{h qq') that we find are in the ballpark of B™'^^{h? qq') 3 x 10~^ 
and B'^'^'^{H^ /A^ ^ gg') ~ 1 x 10~^ for a non-fine-tuned scenario, while B^^^{h^ 
qq') 2 X 10^^ and B^^^{H^ /A'^ ^ gg') ~ 8 x 10^^ for a fine-tuned scenario. From the 
analysis of Ref. [34] we expect that even with these impoverished MSSM rates the number 
of FCNC events of that sort should be non- negligible at the LHC. 

We have already mentioned that our results disagree some orders of magnitude with 
recent estimates presented in the literature [35]. In fact, in our analysis we cannot accom- 
modate a branching ratio at the level of B(h qq') ~ (20 — 30)% for any of the decays 
( |1.1| ), as claimed by these authors. We find such values incompatible with a rigorous MSSM 
analysis of these decays correlated with the branching ratio of b sj. Even though we 
have detected the existence of corners of the MSSM parameter space where a Higgs boson 
FCNC branching ratio can barely reach the 10% level (cf. the narrow windows in Fig. ^ 



^See e.g. Ref. [34] for a detailed analysis of the number of Higgs boson FCNC events produced at the 
LHC in a different situation corresponding to the general 2HDM. 
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and Fig. 0), we insist once more that they should be considered rather unhkely as they are 
associated to fine tuning of the parameters. Moreover, in contrast to these authors, we 
find that it is the hghtest CP-even state, h^, the one that could have the largest FCNC 
branching ratio. Thus, as already advanced in Ref. [34], we believe that the authors of 
Ref. [35] have overestimated by a significant amount the value of B^^^[h qq') for the 
three Higgs bosons of the MSSM. 

To conclude, we have presented a first realistic estimate of the branching ratios of 
the Higgs boson FCNC decays within the MSSM, assuming that the SUSY-QCD 

corrections can be as large as permitted by the experimental constraints on B{b 57). 
We have carried out a systematic and self-consistent maximization of the branching ratios 
( |1.2| ) taking into account this crucial experimental constraint. At the end of the day the 
results that we obtain, especially for the lightest CP-even Higgs boson of the MSSM, are 
fairly large: B'^^''{h^ ^ qq') ^ 10"^ - 10"^. These MSSM rates turn out to be between 
three to four orders of magnitude larger than the maximum SM rate (|2.1| ), but not five or six 
orders as naive expectations indicated. Whether this branching ratio is measurable at the 
LHC [8,9] or at a high energy e+e~ Linear Collider [10] can only be established by means 
of specific experimental analyses. However, on the basis of related studies in the general 
2HDM [34] and from ongoing work in the MSSM [41], we can foresee that an important 
number of FCNC events ( |1.1D can be potentially collected at the LHC. They could play 
a complementary, if not decisive, role in the identification of low-energy Supersymmetry. 
In this paper we have dealt only with the maximum rates induced by the SUSY-QCD 
sector of the model. A more detailed analysis - including the SUSY-EW sector and the 
computation of the aforementioned production rates - will be presented in a forthcoming 
publication [41]. 
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